The normal phase of an imbalanced Fermi gas 
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Recent experiments on imbalanced Fermi gases have raised interest in the physics of an impurity 
immersed in a Fermi sea, the so-called Fermi polaron. In this letter, a simple theory is devised to 
describe dilute Fermi-polaron ensembles corresponding to the normal phase of an imbalanced Fermi 
gas. An exact formula is obtained for the dominant interaction between polarons, expressed solely 
in terms of a single polaron parameter. The physics of this interaction is identified as a signature 
of the Pauli exclusion principle. 
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Quasi-particles are generic emergent properties of 
many-body systems that simplify the description of com- 
plex interacting ensembles of particles. This concept is 
probably one of the most important in quantum physics 
since it lies at the foundation of fields as diverse as chem- 
istry of dilute solutions, dressed atom theory in atomic 
physics or band theory in solid state physics. Recently, 
experiments on spin imbalanced ultra-cold Fermi gases 
[ll-y] have highlighted once again its importance by show- 
ing that the main features of the phase diagram of these 
systems could be understood quantitatively from the 
properties of an impurity immersed in a Fermi sea of spin 
polarized atoms, the Fermi polaron 043- 1^ shown in 
particular that the quasi-particle arising from the inter- 
action between the impurity and the surrounding Fermi 
gas could be described with great accuracy by assum- 
ing that a single particle-hole pair is excited 0, 0| . The 
single-particle properties of the Fermi-polaron have been 
characterized experimentally and theoretically and are 
now well understood. For instance, at unitarity, where 
the scattering length between the majority and minority 
spins is infinite, the chemical potential of the impurity 
is shifted by Hp = Afii where A = —0.61 and /ii is the 
chemical potential of the majority 0, [^, H, . Similarly, 
the effective mass is found to be close to the bare mass 
m, with m* = 1.20to for recent experiment 0, [III [l^ . 
close to the theoretical values obtained from variational 
or Monte-Carlo calculations 0, H, B 

More generally, the Fermi-polaron is a good descrip- 
tion of an impurity immersed in a Fermi sea around uni- 
tarity and for "attractive" (a < 0) interactions where 
A and m* have also been calculated with great accuracy 
0, H, ■ Interestingly, the dressed impurity undergoes a 
transition from a fermionic polaron to a bosonic molecule 
at l/kpia ~ 0.9 H, [13, where kpi = {Gn^ni) 

is the Fermi wave-vector of the majority species gas of 
density ni. This transition reflects into the collective be- 
havior of an ensemble of impurities. In particular, in the 
fermionic sector l/kpia < 0.9, pioneering Fixed Node 



Monte-Carlo simulations have shown that, for a small 
concentration of minority fermions, the equation of state 
of an imbalanced normal Fermi gas with two spin species 
noted a = 1,2 and densities Ua- could be fitted by a 
Landau-Pomeranchuk law 
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where x = n2/rii, EpQi is the energy of a single- 
component (majority) Fermi gas with density ni and F 
describes interactions between polarons 0, 
and F are functions of 1/kpia. This Fermi liquid pic- 
ture is supported by the absence of vortices in rotation 
experiments indicating a normal state [l|. By contrast, 
it was noted recently that experimental data could be 
fitted with great accuracy by a grand-canonical equation 
of state 
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which apparently describes a mixture of two ideal Fermi 
gas of polarons and majority atoms However, the 

presence of a /ii dependence of in the polaron part of 
the equation of state implies a coupling between the two 
gases, and the two equations of state can be reconciled by 
noting that expressed in the canonical ensemble, Eq. ([2]) 
indeed yields Eq. ^ with F = 5A^/9 ~ 0.2 at unitarity, 
close to the Monte Carlo value F 0.14 

In this letter, we show that the equation of state of 
the normal phase follows the phenomenological expan- 
sion ll]). Moreover, we argue that the relationship be- 
tween F and A is exact and can be generalized to the 
full BEC-BCS crossover: indeed, we will show that that 
the parameter F is solely a function of the single polaron 
chemical potential and is given by 
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where fip is computed in the low impurity concentration 
hmit where /zi = Epi = h^kpi/2m. Finally, the study of 
the BCS (Bardeen-Cooper-Schrieffer) regime correspond- 
ing to small and negative value of a allows us to clarify 
the origin of the a;^ term in Eq. ([1]). We attribute it to 
a modification of the singie-polaron properties due the 
Pauli blocking created by the presence of the minority 
Fermi sea and overruling density-mediated polaronic in- 
teractions [i^l which contribute to the higher-order x'^^^ . 

The starting point of our demonstration is the cele- 
brated Luttinger sum-rule, stating that if a many-body 
fermionic system can be analytically connected to an 
ideal Fermi gas 27[, then it possesses a Fermi sur- 



face where the momentum distribution is discontinu- 
ous and which encloses a volume depending on density 
only 2l|, |22| . More quantitatively, the Fermi surface is 
given by the wave- vectors kp^ solutions of the equation 
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where ^^^^ = h k /2m — /ig. and S^. is the self-energy of 
spin a particles. By definition, the single polaron chemi- 
cal potential fip is obtained for vanishingly small impurity 
densities n2 , and is thus solution of the equation [Sf 



^ip = E2(w = 0, fc = 0; Hi,H2 = Mp), 



(5) 



and depends on fii only. Since we consider the situation 
of dilute polarons corresponding to a small concentration 
of impurities, the minority Fermi sea remains small and 
fi2 can be expanded in the vicinity of fip. Let us assume 
for the moment that S2 is analytic in fj,2 and kp2, or 
by rotational invariance: Expanding Eq. ([5|) up to 4th 
order we thus get [11] 
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with Sfj,2 — fJ-2 — fJ-p- The equation of state of the dilute 
impurity gas is obtained from the leading order terms, 
i.e. the first three terms in Eq. ([6|), 
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is the usual definition for the effective mass of a quasi- 
particle. Using Luttinger sum-rule we know that kp2 — 
(67r^n2)^/'^ and combined with Gibbs-Duhem relation 
df^.P = iii, Eq. d?]) leads to a pressure P{pLi) identi- 
cal to Eq. To convert this equation of state in the 
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FIG. 1: Diagrammatic representation of the second order per- 
turbation theory. Integration over frequencies allows one to 
decompose the leftmost diagram into two time-ordered dia- 
grams. Since /i2 is negative, inner minority lines traveling 
backward in time are forbidden, and the rightmost term van- 
ishes. 



canonical ensemble, we use the relationship 
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where Ep2 — h'^k'^p2/2m* and we have neglected higher 
order terms in n2 appearing when taking the derivative 
of m* with /ii. Making use of the definition of the grand 
potential —PV = E — jJ-iNi we finally get Landau 
Pomeranchuk law ([IJ with F given by ([3]) . 

We now verify that Eq. ([3]) is not altered when higher 
orders in Eq. ([6]) are included. Indeed, replacing (5/i2 by 
its leading order expression, terms neglected in Eq. ([6]) 
give rise to a k'^2 contribution to Sii2{kp2). From Gibbs- 
Duhem relation this gives rise to a term cx (/i2 — Mp)''^^ 
in Eq. ([2]), hence a x^/^ contribution to the energy. For 
vanishing x, this term is therefore negligible against x'^ 
and does not contribute to the value of F: this argument 
proves that, provided analyticity conditions are fulfilled, 
Eq. ^ gives the exact value of F. 

We now provide evidence for the analyticity of E2. To 
do so, we make use of a time ordered diagrammatic ex- 
pansion of the self-energy illustrated in Fig. ([Ij [23]. For 
each diagram, a line going forward (backward) in time is 
associated with a 9{£_k,a) {d{—£,k,a)), with 9 the Heavi- 
side step function, and contributes to £,k,a- to the energy 
denominator. In addition, the incoming (outgoing) mi- 
nority line contributes to uj (— w). The main point of the 
argument is the negativity of /ip, and of ^2 for small im- 
purity concentration. Indeed, in this case, ^k^2 is always 
positive, which implies that the Heaviside functions asso- 
ciated with impurities traveling backward in time vanish. 
As a consequence, diagrams containing an impurity loop 
do not contribute to the self energy, and similarly, the 
inner part of the 'main' impurity line cannot travel back 
in time. The denominators are therefore always strictly 
positive and this absence of pole guarantees the analyt- 
icity of E2 . This can be interpreted physically by noting 
that the minority Fermi sea would be empty at these neg- 
ative chemical potentials for vanishing interaction. The 
creation of minority fermions is therefore only triggered 
by interaction processes with the majority component. 
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The above ideas are best illustrated by going to the 
BCS weak coupling limit a — J> 0~, where exact perturba- 
tive calculations can be performed. The gas of fermions 
with two spin-species is described by the Hamiltonian 



fc,(7 
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where Sk = % k"^ /2m, V is a quantization volume, and 
Cfe,(T annihilates a fermion of spin a and momentum fc. 
The zero-range interaction potential in Eq. (jlip suffers 
from ultraviolet divergences that are cured by imposing a 
cutoff kc in momentum space. The Lippmann-Schwinger 
formula then relates the bare coupling constant g to the 
scattering length, 
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Building on the Luttinger equation ^ relating ^2 and 
kF2 for the minority fermions, we wish to determine the 
equation of state P{iJii) in the strongly imbalanced case 
with /xi > and ^2 < 0. The self-energy S2 is calculated 
perturbatively in powers of g. In addition, Eq. (jl2p is 
used to expand the resulting expressions again in powers 
of a. The renormalizability of the model (jlip imposes 
that ultraviolet divergences cancel out for each order in 
a, and the cutoff kc is eventually taken to infinity. 

The first order is given by the usual Hartree diagram, 
T,[^\io,q) = (5/67r2)(2m^i/;i^)3/2. We write the second 
order using the time ordered diagrams displayed in Fig.[TJ 
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where the minority travels partially backward in time in 
the second term and always forward in the first one. As 
stated earlier, the negative minority chemical potential 
implies that Cq+q'-fe,2 is positive and the second term 
of Eq. (jl3p thus vanishes in accordance with our general 
rule that backward travel is suppressed. Moreover, for 
the remaining first term in Eq. ()13p . the denominator 
does not vanish as long as w < — /i2, and the self-energy 
can be freely expanded with respect to [I2 and q at a; = 0. 

Using the complete self-energy I]2^-'-J-v(2) 
to calculate with the result 
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Using Eq. ([3]), we see that up to 3rd order included, 
the interaction parameter F should read 
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It is illuminating to check the weak coupling predic- 
tion ((T5|) for the interaction by a direct calculation of 
the ground state energy using the standard Rayleigh- 
Schrodinger perturbation theory. We first discuss the 
energy of a single polaron £'poi(<7). The unperturbed 
state is then an impurity with momentum q immersed 
in a Fermi sea of majority atoms. The first order cor- 
rection to the energy is the mean- field correction gni, 
while the next order correction involves the excitation of 
particle-hole pairs out of the majority Fermi sea. By def- 
inition of Epai(q) the energy of the system is given by 
E = Efgi + Epoi{q) with 
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where the majority momenta q' and fc' satisfy the con- 
ditions q < kpi (*) and k > kpi (ii) imposed by Pauli 
exclusion principle. 

We switch now to an ensemble of impurities, in which 
case two ideal Fermi gases with Fermi wavevectors kpi 
and kF2 constitute the unperturbed ground state with 
energy EpG,i + Efg.2- The energy takes the form 
E{ni,n2) - EpG.i+E, 
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with the previous restrictions (i), (ii), complemented by 
q < kF2 {in), and \q + q' — k'\ > kF2 {iv), where the last 
two conditions are imposed by the Pauli exclusion prin- 
ciple in the presence of the minority Fermi seas. Except 
for the constraint (iv) E would simply be J2q<kF2 ^poi('3') 
which constitute the energy of an ideal gas of polarons 
with a dispersion relation Epoi{q). However, we can re- 
cover this term explicitly by expressing (iv) in terms of 
its complementary domain (v) |qr-|-q' — fc| < kF2, in which 
case we can recast Eq. (fTT)) as 
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The first term in Eq. corresponds to an ideal gas of 
polarons and contributes to the x and x^^^ scaling terms 
in Eq. ([1]) , that is to A and to* . The second term describes 
the effect of Pauli blocking due to the minority Fermi sea 
on the formation of the polaron. A careful analysis of its 
behavior for low kF2 shows that it scales as x"^ and thus 
gives the effective interaction F between polarons. 

The complete calculation of third order corrections is 
lengthy but straightforward. In the limit a: <C 1, one 
finds again Eq. ([T]) for the ground state energy together 
with an interaction parameter F arising again from Pauli 
blocking and identical to Eq. p^ . 



4 



0.4 



0.3 



n. 0.2 



0.1 




0.0 



FIG. 2: Variations of F in the crossover and comparison with 
Monte-Carlo simulations. Solid: Eq. ((Sjl with calculated 
variationally using @|; Dotted: Third order diagrammatic ex- 
pansion (Eq. (|15|l ): Circles: Mont-Carlo simulations fl^ . 



by investigating higher orders or by making use of non- 
perturbative schemes. In fact, assuming further analyt- 
icity, the low density expansion performed here can in 
principle be extended to any order in x. The coefficients 
of the expansion are then expressed solely in terms of the 
single-polaron self-energy. Other open questions include 
the extension of our results to the one-dimensional situa- 
13, 25 1 and to the case of repulsive interactions [1^ . 
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The argument presented above makes a strong case 
for a interaction between polarons. However, notic- 
ing that s-wave interactions gives a scaling and p- 
wave a subleading x^/^, this may seem to contradict the 
fermionic nature of polarons. On the other hand, Fermi 
liquid theory does not forbid alike particles to interact, 
and the corresponding interaction is in fact not necessar- 
ily short-ranged. This paradox can be solved by noting 
that polarons have fermionic statistics at large distances 
and are composite objects at shorter distances. From this 
structure, they acquire an internal energy ^.p = A/ii. 
This single polaron energy is held fixed in the grand- 
canonical ensemble and is not modified by the presence 
of other impurities. By contrast the internal energy de- 
pends on the minority concentration in the canonical en- 
semble through Pauli blocking which yields the inter- 
action in Eq. Based on these arguments, it is prob- 
ably not surprising to find that F is solely a function of 
the internal energy as given by Eq. ([3]). 

Finally, in Fig. [2l we compare our prediction Eq. (|3]) 
where /Zp is calculated using the variational scheme pre- 
sented in [il with the third order expansion Eq. (fT5|) as 
well as Monte-Garlo data [1^]. As expected, we observe 
that the perturbative expansion and the non perturba- 
tive result coincide for a — ^ 0~. In the strongly inter- 
acting limit we observe that our result follows the same 
trend as the Monte-Garlo simulation, with in particular 
the presence of a maximum of F close to l/kpia ^ 0.5. 

In conclusion, we have demonstrated that in the low 
impurity concentration, the canonical equation of state of 
a spin imbalanced system could be described by a Lan- 
dau Pomeranchuk energy. Quite surprisingly, we have 
shown that the interaction parameter F was related to 
single impurity properties. Several extensions of this let- 
ter are worth exploring. From experimental data, it ap- 
pears that Eq. ([2]) is valid on a wide range of impurity 
concentrations (up to x = 0.5 at unitarity). This sur- 
prisingly large validity domain remains to be understood 
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